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Abstract — In this paper, a closed-form solution minimizing the 
Godard or Constant Modulus (CM) cost function under the 
practical conditions of finite SNR and finite equalizer length 
is derived. While previous work has been reported by Zeng 
et aL, IEEE Trans. Information Theory. 1998, to establish the 
link between the constant modulus and Wiener receivers, we 
show that under the Gaussian approximation of intersymbol 
interference at the output of the equalizer, the CM finite-length 
receiver is equivalent to the nonblind MMSE equalizer up to 
a complex gain factor. Some simulation results are provided to 
support the Gaussian approximation assumption. 

I. Introduction 

The transmission of information over frequency selective 
digital communications channels is subject to intersymbol 
interference (ISI). Equalization has proven to be an effective 
means for removing ISI. The commonly used approach em- 
ploys a known sequence of training symbols and the equalizer 
coefficients are then adapted by using some adaptive algorithm 
so that the output of the equalizer closely matches the training 
sequences by minimizing the mean-squares-error (MMSE) 
criterion. However, in many high data rate bandlimited dig- 
ital communications systems, the transmission of a training 
sequence results in a significant reduction in the effective 
communication link data rate. Therefore, instead of using a 
training sequence, only some statistical or structural properties 
of the transmitted and received signal can be exploited in a 
blind process to adapt the equalizer i.e. blind equalization. 
One of the simplest and the most effective blind equalization 
schemes is the constant modulus algorithm (CM A) [1], [2], 
which uses constant modularity as the desired property of 
the output of the receiver. A great deal of research about 
convergence behavior has been reported but a lack of com- 
prehension about the global convergence of these algorithms 
has limited their utilization. Global convergence has been 
proven under ideal conditions [1], [2], [3], [4] (e.g., noise-free 
and doubly infinite equalizer) providing the inverse transfer 
function of the channel. However, this assumption can be 
excessive and the channel noise must be considered in the 
design of the equalizer. Recently, Zeng et al. [5] have proposed 
interesting results on the relashionships between blind and 
Wiener receivers in the noisy case. They observe that if a 
Wiener equalizer reaches an acceptable MSE performance, 
there exists a CM equalizer in its immediate neighborhood 
which is approximately a scaled version of the Wiener receiver. 
Nevertheless, never has a closed-form solution for the CM 



equalizer been proposed in the noisy case. In this paper, a 
closed-form solution for the finite-length CM equalizer for a 
non-minimum phase channel in the presence of additive white 
Gaussian noise is given. The CM equalizer is shown to be 
equivalent to the MMSE equalizer up to a complex gain factor. 
In order to illustrate the theoretical result, a simulation for 
4PSK signals is performed, where a comparison of MMSE 
and CM equalizers, using gradient stochastic algorithms, is 
implemented. 

II. A MODIFIED EXPRESSION OF THE CM COST FUNCTION 

We consider a discrete-time baseband transmission model. 
Complex symbols a n are independent, identically distributed 
(i.i.d.) random variables with zero-mean, variance a\ and sub- 
Gaussian distribution (K(a) < 0). Symbols are sent through 
a discrete-time equivalent channel with coefficients hi. The 
output of the channel is corrupted by ISI and can be expressed 
by 

L 

i=o 

where b n is a zero-mean, white Gaussian complex and 
circularly-symmetric random process with variance a 2 and 
L is the channel memory. Considering a block of M output 
symbols and assuming that the channel is time-invariant over 
this block, Q can be written more compactly as follows 

Y n = HA n + B n (2) 

where we define the column vector 

A [ ~l T a tr>M+L—l 

A n — [ a n a„_i . . . a„_( M+£ _ 2 ) J t 

of the input symbol sequence and the additive noise 

vector B n = [ b n b n ^ ... 6„_ (A /-i) ] T G C M . 

H is the M x (M + L — 1) channel Toeplitz matrix. 

Denote superscript T as the transpose operator. If 

C = [ c ci ... cm-i ] T e C M is the tap-weight 

vector of the equalizer of length M, the output of the 

equalizer z n can be represented by 

z n = C H Y n (3) 

where superscript H is the Hermitian transpose operator. The 
CM cost function [1] with index p = 2 is defined by 

= E{(\z n \ 2 - R 2 f} (4) 



where R 2 is a real positive dispersion factor. A necessary 
condition for the minimization of @ is given by 

V C D^ = E{Y n Y^C(\z n \ 2 - R 2 )} = (5) 

where we define the gradient of at C as the row vector 
Vc-t^ 2 - 1 . Then, one can write the following relation 

C H V C D^ = E{\z n \ 2 (]z n \ 2 - R 2 )} = (6) 

Finally, the condition necessary to reach a minimum of the 
CM cost function involves the equality between the second 
and fourth order moments of z n up to a positive gain R 2 

E{\z n \ 4 } = R 2 E{\z n \ 2 } (7) 

Considering now the development of (|4}, we have 

= R 2 (R 2 - E{\z n \ 2 }) + E{\z n \ 4 } - R 2 E{\z n \ 2 } (8) 

Then, a cost function equivalent to subject to the neces- 
sary condition (0 can be defined by 

= R 2 (R 2 - E{\z n \ 2 }) (9) 

Thus, the minimization of under the constraint @ is 
equivalent to the minimization of the CM cost function D^ 2 \ 

III. Finite-length equalizer minimizing the CM 

COST FUNCTION 

From (0 and 0, the output of the equalizer is given by 

z n = S H A n + C H B n (10) 

where S denotes the column vector that contains the coef- 
ficients of the combined impulse response sj of the channel 
and equalizer as illustrated in Fig. ^ Let w n and x n denote 

b„ I I w„ 
°l 1 
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combined impulse response of the channel and equalizer 

Fig. 1. Channel- equalizer block diagram. The cascade channel-equalizer is 
known as the global response. 

the noise and the residual interference after filtering by the 
equalizer, respectively. In this paper, our approach differs from 
work on channel deconvolution because we assume that the 
residual interference is not null at the equalizer output, x n ^ 0. 
So, the output of the equalizer can be written 

Z n SyQj n — v -\- X n -\~ W n (11) 

where the receiver gain s v is defined as C H He„, and e„ 
denotes the vth unit row vector (contains only one non zero 
for a particular value of delay v). The Kurtosis of z n can be 
expressed by 

K{z) = E{\z n \ 4 } - 2(E{\z n \ 2 }) 2 - \E{zl}\ 2 (12) 



where the last term is null when z n is circular, e.g. when a n is 
a complex-value so that ^{a 2 } = 0, leading to E{z 2 } = 0. 
Since symbols a n are assumed Lid. and independent of the 
noise components, the Kurtosis of z n can be written as a sum 
of separate Kurtoses 

K(z) = |s„| 4 /v (a) + K{x) + K(w) (13) 

When the number of taps of the equalizer is sufficiently 
large, the pdf of the residual interference x n is reasonably 
well modeled by a Gaussian distribution. Although the ISI 
distribution for a discrete input can never be Gaussian, it is 
commonly accepted that the ISI at the output of an MMSE 
equalizer may be estimate by a Gaussian distribution [6]. 
This assumption is based on the central limit theorem and is 
discussed in section Hvl Thus, since the Kurtosis of a Gaussian 
random variable is null, we have 

K(z) = \s„\ 4 K(a) (14) 

Substituting (Till in il2\ gives 

|,s„| 4 ^( a ) - E{\z n \ 4 } - 2(E{\z n \ 2 }) 2 (15) 

So, using necessary condition yields 

\s„\ 4 K(a) = R 2 E{\z n \ 2 } - 2(E{\z n \ 2 }) 2 (16) 

The variance of the receiver output z n can be expressed 
according to Q as 

E{\z n \ 2 } = C H R yy C (17) 

where R yy is the output autocorrelation matrice of the received 
data sequence defined by 

R yy = E{y n y*J = a 2 b I M xM + o 2 a BH H (18) 

Denote * superscript as the complex conjugate operator. In 
order to obtain a closed-form for the equalizer C, using the 
gain s u we derived the Schwarz inequality which allows us to 
write 

| s „| 2 = \C H He„\ 2 < {C H R yy C) (e^H H R^He,) (19) 

with equality if and only if 

ae" H H = C H R yy (20) 

where a is a complex coefficient. The above inequality follows 
from 1 (assuming that R yy is invertible). Let us define 

Q v = H H R^ He y (21) 

According to Mil , expression dl9l can be rewritten by 

|s,| 2 < n u E{\z n \ 2 } (22) 

Considering the sub-Gaussian case (K(a) < 0), K(a) = 
— \K(a)\. Using d!6i and i22\ . we show that the necessary 

'we consider the unique factorization of the positive-definite matrix R vy 
in the form R yy = G 1 / 2 G H / 2 and then the decomposition C H He v = 



condition Q to obtain a minimum of the CM cost function 
involves 

i?2 



E{\z n \ 2 }< 



2-\K(a)\ni 



(23) 



Substituting now (I23> in l|9}, we can see that is lower 
bounded by 



D {2) > (i? 2 ) 2 (l 



-) 



(24) 



2-\K(aW u 

So we deduce that will be minimum if and only if 

C H = ae^ H H Ryy (25) 

Since the minimization of with respect to the necessary 
condition |7} is equivalent to the minimization of (section 
II), the transfer function i25\ corresponds to a global mini- 
mum. Then requiring the equality of (1231 . it results that the 
output power is bounded by R2/2 < E{\z n \ 2 } < R2 which is 
consistent with [5] in the complex case. To find a, we replace 
C by its value d25l in dl7l and use (I23> . Then, the solution is 
given by 

12 1 ^2 



0,2- \K{a)\9l 



(26) 



Finally, the finite-length equalizer minimizing the CM criteria 
is defined up to a phase ambiguity by 

C H = \a\e* H H (a 2 b I MxM +a 2 a HH H )- 1 exp(j0), 6 € [0; 2tt] 

(27) 

Hence recalling the finite-length MMSE equalizer minimizing 
the mean square error which is a causal Wiener filter given by 



H 



e^H H a 2 a {aiI MxM + KHH 



H\-l 



(28) 



And using (I27> and d28i . one can easily obtain a relation 
between the MMSE and CM equalizers 



a 



H = exp{-jO) ^ 0e[o . 2lr] (29) 
M erf 



Since MMSE and CM equalizers are identical up to a complex 
factor, it is theoretically possible to reach the MMSE equalizer 
performance with a blind receiver composed of a CM equalizer 
followed by a correction gain and a phase rotator, as depicted 
in Fig. 
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equivalent to the MMSE equalizer 
Fig. 2. CM equalizer equivalent to the MMSE equalizer. 



IV. Statistical measure of Gaussian distribution 

The output of the blind and MMSE receiver is shown to be 
theoretically the same. Consequently, the residual interference 
at the output of the two receivers will be similar and we 
can equivalently study the Gaussian nature of the residual 
interference at the output of the optimum finite length MMSE 
equalizer. A commonly used approach [7] is to plot the data 
against a theoretical normal distribution in such a way that if 
the data set exhibits the properties of a normal distribution 
the points should lie in a straight line. Deviations reflect 
miss-matches between the data distribution measured and 
normal distribution. Specifically if < x^ 2 \.. < x^ N > 
denote the N ordered observations x n , n = 1,...,N of the 
residual interference in (II U and q( l=1 <••■ ' N > the corresponding 
quantiles, mean the fraction of points below a given value, of 
the standard normal distribution [8], the N points 
i = l,...,N define the normal probability plot. The dis- 
tribution of the ISI output is also reported and compared 
to a theoretical gaussian pdf. We have considered various 
numbers of equalizer coefficients and results are given by 
generating N=250000 symbols for two signal to noise ratios 
(SNR) defined by SNR = al\\h\\ 2 / a% . The channel is time- 
invariant and given by the following set of real coefficients 
[0.5679 -0.1136 0.5849 0.1124 0.556]. The distributions are 
normalized in location and scale. The normal probability plot 
in Fig. |5J-Left shows a strongly linear pattern in the center of 
the data. The quality of the fit is evaluated through the Pearson 
correlation coefficient defined by 



1 N 



(30) 



which measures the correlation between the two distributions, 
a value of +1 indicating a perfect positive linear relationship. 
We observe that the first and the last few points deviate 
from the reference line since the number of tap gains M of 
the equalizer is not large enough. It is noticeable in Fig. |3j- 
Right-a with distortions of the empirical pdf relative to the 
normal distribution. The same observations are related in 
Fig. |4}Left /|4}Right for an SNR of 15dB. However, in all 
cases, the tendency towards normality is obvious. A further 
characterization of the distribution relative to the standard 
normal one includes the kurtosis measure defined as 



1 N 

n=l 




(31) 



As indicated in Fig. [3] and Fig. |4] the kurtosis is still near 
zero which is the key assumption of the task. We can quite 
conclude that for a reasonable SNR/equalizer length, relation 
H4i is achieved. 

V. Simulation results 

To confirm our theoretical results, we propose to compare 
the response of the proposed blind receiver as depicted in 
Fig. 13 with that of the Wiener solution. If we assume that 
emitted symbols are known at the receiver, we can use a 
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Fig. 3. Gaussianity measure. Left: Normal probability plot. The points should 
fit the reference line to check that the data set is normally distributed. Right: 
Probability density function. We observe the similarity between the empirical 
pdf of the data and the theoretical Gaussian pdf. SNR=20dB a) M=l 1 b) 
M=21 c) M=41. 
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Fig. 4. Gaussianity measure. Left: Normal probability plot. The points should 
fit the reference line to check that the data set is normally distributed. Right: 
Probability density function. We observe the similarity between the empirical 
pdf of the data and the theoretical Gaussian pdf. SNR=15dB a) M=ll b) 
M=21 c) M=41. 



data-aided least mean square (LMS) algorithm that converges 
towards the Wiener solution. The tap weight update vector of 
the equalizer is as follows 



CUs^) = CUs^ - 1) - v{z n - a n )Y* 



(32) 



where /i is a step-size parameter. Using a stochastic gradient 
of the CM cost function (|4} with respect to the tap weight 
vector of the blind receiver, each iteration of the algorithm 
involves updating the following relation 



C(n) = C(n - 1) - fJ-z„(\z n 



R2)y: 



(33) 



From relation ( I29> . we need an estimation of \a\, denoted |d|. 
This value can be obtained blindly from an estimate of the 
power at the equalizer output. Considering J23t when the CM 



cost function is minimum 

E{\Zn\ 2 } = 



i?2 



2- \K(a)\n 2 
from (I26> and J34> . we can write 



(34) 



(35) 



with g\ being estimated recursively from the following rela- 
tion 




a 2 z {n) = Xa 2 z (n - 1) + (1 - \)\z n f 



(36) 



Finally, the phase rotator can be realized from a second-order 
carrier-phase tracking loop operating in a decision-directed 
mode. Simulations consider a 4PSK baseband model where 



1, K(a) 



1 and i?2 



1. The channel is the 



same set of real coefficients as that of section II VI We have 
considered a high (M — 41) number of complex tap gains 
for the equalizers in order to avoid unstable phenomena due 
to finite-length equalization constraints [9]. The center-tap is 
initialized to 1 and the step-size /j, is set to 0.001. Results are 
given after the transmission of 50000 symbols (the algorithm 
performs one update iteration per symbol). As shown in Fig. EI 




Finally we would like to emphasize that the proposed solution 
clearly deserves further investigation. 
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Fig. 5. Comparison of data-aided (LMS) and blind equalizer coefficients after 
parameter convergence for SNR=20dB and SNR=15dB. The blind receiver 
response is very close with that of the MMSE equalizer. 



the response of the blind receiver including the CM equalizer 
and the complex gain compensation is very close to that of 
the LMS equalizer for both SNR. Consequently, the bit error 
rate performance of the two receivers will be similar. It should 
be noted that even if asymptotic performance is the same, the 
data-aided LMS algorithm exhibits faster convergence than the 
blind algorithm. 



VI. Conclusion 

In the present paper a closed-form solution for a finite length 
CM equalizer under a non-minimum phase channel in the 
presence of additive white Gaussian noise is presented. The 
analytical resolution is based on the Gaussian approximation 
assumption of intersymbol interference at the output of the 
blind receiver. Statistical measures are used to validate the hy- 
pothesis. In the case of a reasonable length equalizer and finite 
SNR, the empirical output distribution tends to the normal one 
and the zero kurtosis assumption is valid. Under this condition, 
the CM equalizer is shown to be equivalent to the MMSE up 
to a complex gain factor, which is consistent with the relation 
between the CM and Wiener receivers proposed by Zeng et 
al. [5]. Therefore the CM receiver can theoretically reach the 
same asymptotic performance as the nonblind receiver. The 
simulation results herein corroborate the theoretical analysis. 



